An ab initio theory of double odd-even mass differences in nuclei 
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Abstract. Two aspects of the problem of evaluating double odd-even mass differences D 2 in semi-magic nuclei 
are studied related to existence of two components with different properties, a superfluid nuclear subsystem and 
a non-superfluid one. For the superfluid subsystem, the difference Di is approximately equal to 2A, the gap A 
being the solution of the gap equation. For the non-superfluid subsystem, D 2 is found by solving the equation for 
two-particle Green function for normal systems. Both equations under consideration contain the same effective 
pairing interaction. For the latter, the semi-microscopic model is used in which the main term calculated from 
the first principles is supplemented with a small phenomenological addendum containing one phenomenological 
parameter supposed to be universal for all medium and heavy atomic nuclei. 
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1 Introduction 

Recently, a progress has been made in the 'ab initio ' the- 
ory of nuclear pairing, first, by the Milan group I11J2J31 . a 
bit later by Duguet et al. [4 5 1, and finally by the Moscow- 
Catania group [6789 10l . By 'ab initio' we mean a the- 
ory of Brueckner type starting from a free AW-potential, 
not a real ab initio approach based on the QCD theory. 
The quotes indicate the fact that, in any of the cited refer- 
ences, the authors deal with a calculation which is not com- 
pletely ab initio even at such lowest level. Indeed, the self- 
consistent mean field used in each of them is calculated by 
means of an Energy Density Functional (EDF) containing 
phenomenological parameters. Solving the BCS equation 
for the gap A with a realistic AW-potential for the pair- 
ing interaction is the main ingredient of all of them. This 
is exactly the prescription of the simplest version of the 
Brueckner theory as the ladder diagram summation is al- 
ready performed in the gap equation itself. 

In the first paper of the Milan series, the BCS gap equa- 
tion for neutrons with the Argonne v 14 potential was solved 
for the nucleus 120 Sn, the Saxon- Woods Shell-Model ba- 
sis with the bare neutron mass m* = m. Rather optimistic 
result was obtained for the gap value, zIbcs = 2.2 MeV 
which is bigger, but not dramatically, than the experimen- 
tal one, zf exp ^1.3 MeV, leaving some hope of achieving 
a good agreement by finding corrections to the scheme. In 
Refs. I2I3L the effective mass m* + m was introduced into 
the gap equation within the Skyrme-Hartree-Fock (SHF) 
method with the Sly4 force 1111 . In this case the effec- 
tive mass m*(r) is space dependent and essentially differ- 
ent from the bare one m. E.g., in nuclear matter the Sly4 
effective mass is equal to m* = 0.7m. In the weak cou- 
pling limit of the BCS theory, the gap is exponentially 
dependent, i.e. A cc exp(l/g), on the inverse dimension- 
less pairing strength g - ni"V e $k? I n 2 , where < V e fs is the 
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effective pairing interaction. Therefore, a strong suppres- 
sion of the gap takes place in the case of m* < m. The 
value of zIbcs = 0.7 MeV was obtained in Ref. [2] and 
^bcs = 104 MeV, in Ref. 0. In both cases, the too small 
value of the gap was explained by invoking various many- 
body corrections to the BCS approximation. The main cor- 
rection is due to the exchange of low-lying surface vibra- 
tions ("phonons"), contributing to the gap around 0.7 MeV 
(2), so that the sum turns out to be A = 1.4 MeV very 
close to the experimental value. In Ref. [3|, the contribu- 
tion of the induced interaction caused by exchange of the 
high-lying in-volume excitations was added either, the sum 
again is equal to A 1.4 MeV. Thus, the calculations of 
Refs. [2 3 1 showed that the effects of m* + m and of many- 
body corrections to the BCS theory are necessary to ex- 
plain the difference of (zIbcs _ ^exp)- In addition, they have 
different sign and partially compensate each other. Unfor- 
tunately, both effects contain large uncertainties. This point 
was discussed in Refs. B7I9I101 . 

To avoid such uncertainties, a semi-microscopic model 
for nuclear pairing was suggested in Refs. 118191101 . It starts 
from the ab initio BCS gap equation with the Argonne 
force Vis treated with the two-step method. The complete 
Hilbert space S of the problem was split into the model 
subspace So and the complementary one S'. The gap equa- 
tion is solved in the model space with the effective inter- 
action "Veff which is found in the complementary subspace 
in terms of the initial AW-potential *V with m* = m. This 
ab-initio term of 'Y^ was supplemented by a small ad- 
dendum proportional to the phenomenological parameter 
y that should hopefully embody all corrections to this sim- 
plest BCS scheme. This parameter is supposed to be the 
same for all heavy and medium nuclei, either for neutrons 
and protons. 

The "experimental" gap value A exp for semi-magic nu- 
clei is usually identified with one half of one of that ob- 
tained from the following double odd-even mass differ- 
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ences: 

D+ n (N, Z) = M(N + 2, Z) + M(N, Z) - 2M{N + 1 , Z), ( 1 ) 

D 2n (N, Z) = -M(N - 2, Z) - M(N, Z) + 2M(N- 1 , Z), (2) 

Z)^(AT, Z) = M(N, Z + 2) + M(JV, Z) - 2M{N, Z + 1), (3) 

D 2p (N,Z) = -M(N,Z-2)-M(N,Z) + 2M(N,Z-l), (4) 

where Z value is magic in Eqs. ©,© and N is magic in 
Eqs. ©,©. More explicitly, in Refs. 1191 101 the experimen- 
tal gap values were defined as follows: 

2AT V (N, Z) = 1 l2{DlfN, Z) + D 2n (N, Z)), (5) 

2A e p xp (N, Z) = 1 /2(D* p (N, Z) + D 2p (N, Z)). (6) 

The accuracy of such a prescription was estimated in Q as 
^0.1 -r 0.2 MeV. Approximately the same accuracy holds 
for the "developed pairing" approximation in the gap equa- 
tion, with conservation of the particle number only on av- 
erage [ 12|, used in [1-10]. 

There is one more physical quantity in semi-magic nu- 
clei which can be evaluated in terms of the same effective 
interaction as the pairing gap. This is the set of the same 
double odd-even mass differences dTJ — (EJ>, but for the non- 
superfiuid subsystem, so that now N is magic and Z ar- 
bitrary in Eqs. ©,© and vice versa in Eqs. ©,©. In 
non-superfluid nuclei, the mass differences, Eqs. ©,©, 
coincide with poles in the total energy E plane of the two- 
particle Green function K( 1 , 2, 3 , 4) for normal systems [ 1 3 1 
in the nn-channel, and Eqs. ©,©, in the /?/?-channel. The 
equation for K in the channel S = 0, L — could be ex- 
pressed in terms of the same effective interaction T^ff as 
the pairing gap. This point was marked in the old paper 
fl4"l . where these differences for double-magic nuclei were 
analyzed within the theory of finite Fermi systems (TFFS) 
1 1 3 1 . In this article, the density dependent effective pairing 
interaction was introduced for the first time and arguments 
were found in favor of of the surface dominance in this 
interaction. 

It is worth to stress that this calculation of mass dif- 
ference for the non-superfluid subsystem is a more rigor- 
ous operation than the identification with the double gap 
A in the superfluid one. The first results of such calcula- 
tions with the use of the semi-microscopic model for the 
effective pairing interaction are presented in lfj"5l . 

2 The semi-microscopic model for nuclear 
pairing 

The general many-body form of the equation for the pair- 
ing gap is ED 

A T = U T G T G S T A T , (7) 

where r = (n,p) is the isotopic index, tl T is the NN- 
interaction block irreducible in the two-particle T-channel, 
and G T (G S T ) is the one-particle Green function without 
(with) pairing. A symbolic multiplication denotes the inte- 
gration over energy and intermediate coordinates and sum- 
mation over spin variables as well. The BCS approxima- 
tion in Eq. © means, first, the change of the block H of 
irreducible interaction diagrams with the free AW-potential 
*V and, second, the use of simple quasi-particle Green func- 
tions G and G s , i.e. those without phonon corrections and 



so on. In this case, Eq. © iss greatly simplified and can be 
reduced to the form typical of the Bogolyubov method, 

A T = -V T x T , (8) 

where 

/ds 
^-G T G S T A T (9) 

is the anomalous density matrix which can be expressed 
explicitly in terms of the Bogolyubov functions u and u, 

^(ri,r 2 ) = ^ M J(r 1 ) I ;J(r 2 ). (10) 

i 

Summation in Eq. (TTOb scans the complete set of Bogolyubov 
functions with eigen-energies E, > 0. 

Then we split the complete Hilbert space of the pair- 
ing problem, S — So + S', where the model subspace So 
includes the single-particle states with energies less than 
a separation energy Eo. The gap equation is solved in the 
model space: 

A T = ^G T GlA T \ So , (11) 

with the effective pairing interaction < V BC ^! instead of the 
block 'V in the BCS version of the original gap equation 
©. It obeys the Bethe-Goldstone type equation in the sub- 
sidiary space, 

«VB« _ ryr + ryTQ Q ryBCS {u) 

In this equation, the pairing effects can be neglected pro- 
vided the model space is sufficiently large, Eo » A. That is 
why we replaced the Green function G S T for the superfluid 
system with its counterpart G T for the normal system. To 
solve Eq. ( fT2l in non-homogeneous systems a new form of 
the local approximation, the Local Potential Approxima- 
tion (LPA), was developed by the Moscow-Catania group. 
Originally, it was developed for semi-infinite nuclear mat- 
ter, then applied to the slab of nuclear matter (see review 
article [ 16|) and, finally, to finite nuclei 06171 . It turned out 
that, with a very high accuracy, at each value of the aver- 
age cm. coordinate R = (ri + r 2 + + r^/A, one can use 
in Eq. (fT2l the formulae of the infinite system embedded 
into the constant potential well U = U(R). This signifi- 
cantly simplifies the equation for < V e ff , in comparison with 
the original equation for A . As a result, the subspace S ' can 
be chosen as large as necessary to achieve the convergence. 
The accuracy of LPA depends on the separation energy Eo- 
For finite nuclei, the value of £o=40 MeV guarantees an 
accuracy higher than 0.01 MeV for the gap A. 

To avoid uncertainties that affect the corrections to the 
BCS scheme discussed above, a semi-microscopic model 
was suggested in Refs. B8I9I10L In this model, a small phe- 
nomenological addendum to the effective pairing interac- 
tion is introduced which embodies approximately all these 
corrections. The simplest ansatz for it is 

^eV r i ■ r 2, r 3 , r 4 ) = 1^(ri , r 2 , r 3 , r 4 ) 




Here p(r) is the density of nucleons of the kind under 
consideration, and y T are dimensionless phenomenological 



NSRT12 



> 





2,0 



Fig. 1. Neutron gap in Sn isotopes 
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Fig. 2. Neutron gap in Pb isotopes 



parameters. To avoid any influence of the shell fluctuations 
in the value of p(0), the average central density p(Q) is used 
in the denominator of the additional term. It is averaged 
over the interval of r<2 fm. The first, ab initio, term in the 
r.h.s. of Eq. (JT3J is the solution of Eq. (TTZl i in the frame- 
work of the LPA method described above, with m*-m in 
the subspace S'. 

In Ref.[9], the above equations were solved in the self- 
consistent /i-basis of the EDF by Fayans et al. B17I18II . Two 
sets of the functional were used, the original one DF3 IfTFI 
and its modification DF3-a 1191 . In the latter, the spin-orbit 
and effective tensor terms of the original functional were 
modified. The results for the pairing gap in three chains of 
semi-magic nuclei are displayed in Figs. 1-3. 

In accordance with the recipe of Ref. [3 1, we represent 
the theoretical gap with the "Fermi average" combination 

A F = £ (2,+ DzW 2(2/' +1 )> (14) 



where the summation is carried out over the states A in 
the interval of |e^— //|<3 MeV. The "experimental" gap is 
determined by the symmetric 5-term odd-even mass differ- 
ences d!),©. 

Let us begin from neutron pairing and consider first 
the tin isotopes, Fig. [T] We see that the BCS gap (y=0) 
is approximately 30% greater than the experimental one. 
Switching on the phenomenological addendum with y=0.06 
makes the theoretical gap values closer to experiment. How- 
ever, the predictions of the two versions of the functional 
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Fig. 3. Proton gap for N = 82 isotones 
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used are significantly different, being much better for the 
DF3-a functional. In particular, the A-dependence of the 
experimental gap is reproduced with a pronounced mini- 
mum in the center of the chain. As the analysis in Ref.[9| 
has shown, this strong difference between results for two 
functionals is due to the strong influence on the gap of the 
high j intruder state I/jh/2- Its position depends essentially 
on the spin-orbit parameters and is noticeably different for 
DF3 and DF3-a functionals. 

In the lead chain, see Fig.[2l the overall pattern is quite 
similar. Again the BCS gap is approximately 30% bigger 
of the experimental one and again the inclusion of the phe- 
nomenological term with y=0.06 gives a qualitative agree- 
ment. Now, the difference between the two functionals is 
much smaller. In this case, the agreement is rather perfect 
for the DF3 functional and a little worse for the DF3-a one, 
but also within limits for the accuracy discussed above. 

Let us go to proton pairing, A^=82 chain, see Fig. [3] In 
this case, the Coulomb interaction should be included into 
the pairing effective interaction, 



Y eff 



(15) 



As it is argued in Ref. [9 1, the bare Coulomb potential could 
be used in this equation with high accuracy. The strong 
Coulomb effect in the gap is demonstrated in Fig. [3] It 
is also explained with the exponential dependence of the 
gap on the pairing interaction. It should be mentioned that 
Duguet and co-authors [5] were the first who inserted ex- 
plicitly the Coulomb interaction into the pairing force for 
protons. Only after inclusion of the Coulomb interaction 
into 'Vetf, we can use the same value of y — 0.06 for pro- 
tons and neutrons. As for tin isotopes, the difference be- 
tween DF3 and DF3-a results is rather strong, now in favor 
of the DF3 functional. This effect is again due to different 
positions of the l/zn/2 level, but now for protons. Overall 
agreement with experiment is for protons worse, maybe, 
because of closeness of some nuclei to the region of the 
phase transition to deformed state. 

As it is argued in Introduction, we may expect to reach 
accuracy for the gap of the order of st (0.1 -r 0.2) MeV. 
In practice, the overall disagreement of the theory with the 

data is J(Jth -^exp^-0.13 MeV for the DF3 functional 
and ^0.14 MeV, for the DF3-a one. For the mass differ- 
ences £>2 these numbers should be multiplied by a factor 
of two. 
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3 Double mass differences in 
non-superfluid subsystems 

The Lehmann expansion for the two-particle Green func- 
tion Kin a non-superfluid system reads, in the single-particle 
wave functions |1)=|tii,/i, j\, m\) representation, [ 13]: 



e-e: 



: iy 



(16) 



where E is the total energy in the two-particle channel and 
£j ~ denote the eigen-energies of nuclei with two particles 
and two holes, respectively, added to the original nucleus. 
They are often interpreted as the "pair vibrations" Il20l . In- 
stead of the Green function K, it is convenient to use the 
two-particle interaction amplitude f : 



K - Kq + K^rKo, 



(17) 



where Kq - GG. The amplitude r obeys the following 
equation IfPJl : 

r = <U + <UGGr, (18) 

where U is the same irreducible interaction block as in 
Eq. (0. Again, within the Brueckner theory, the block 1A 
should be replaced with the realistic potential *V which 
does not depend on the energy. Then the integration over 
the relative energy can be readily carried out in Eq. ([18}: 

/de IE \ IE \ 1— n\— n 2 
_ Gl (_ +e ) G2 (__ e )._L_^, (19) 

where e\^ 2 are the single-particle energies and «i,2=(0; 1), 
the corresponding occupation numbers. As the result, we 
obtain: 

r = <V + <VAR (20) 

Simple transformations in Eq. (TT7T > or (l20l . in vicinity 
of the pole E=E*'~, lead to the following equation for the 
eigenfunctions x s '- 



(E s -ex- s 2 ) X \ 2 = (1 - m - n 2 )J^^V 3 i y 34 . 



(21) 



It is different from the Shrodinger equation for two in- 
teracting particles in an external field only for the factor 
(1 — Til — n<£) which reflects the many -body character of the 
problem, in particular, the Pauli principle. The direct solu- 
tion of this equation is complicated by the same reasons as 
of the ab initio BCS gap equation described in Sect. 2. The 
same two-step method is used. The usual renormalization 
of Eq. (EH transforms it into the analogous equation in the 
model space: 



{E s -ei-s 2 )x[ 2 ={l-ni-n z ) i £ A CV e ff)n *34> 

34 



(22) 



where the effective interaction *V e ff coincides with that of 
pairing problem, Eq. ([T2l . provided the same value of the 
separation energy £0 is used. The next step consists in the 
use of ansatz ([T3l to take into account corrections to the 
Brueckner theory with a phenomenological addendum (~ 

y). 

The double mass differences (|T} — (|4]i are identified 
with the two first solutions E*~ of Eq. ( 1221 . correspond- 
ing to the addition of two particles (holes) to the magic 
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two protons 
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core into the state s\ — e 2 = fi + ~ , where the chemical po- 
tentials fj. +, ~ are defined in a usual way as mass differences, 
e.g., ju+ = E B (N, Z + 1) - E B (N, Z). Then, the energy differ- 
ence in the left-hand side of Eq. (l22l is Es'~-2fi + ~ = 
If we are interested only in these solutions, we may rewrite 
Eq. ( 122b as follows: 



E s -2»=a-2ni)(r'(E s ))\\, 



(23) 



where 



(r(E s ))^ = (T eff )i2 + £V eff ) 
56 



E s -ss-ee 



(24) 



The accent on the sum denotes that the two-particle state 
5=6=1 is excluded. In all formulas (l22l - d24"l > the angular 
momenta of two-particle states 1 12), |34) are coupled to the 
total angular momentum 7=0. 

The set of Eqs. (1231. (1241 1 was solved with the functional 
DF3 only. For lead isotopes, the proton differences are 
displayed in Fig. |4] and the differences D 2 , in Fig. [5] and 
for the tin isotopes, analogously, in Fig. |6jand|7] In simi- 
lar way, the neutron double differences for isotones N=&2, 
are displayed in Figs. [8] and [9] Note that the opposite signs 
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of the quantities Dt and D 2 occur owing to the factor of 
(1 - 2«i) in Eq. d23l >: particles attract each other whereas 
holes, repel. The results of the ab initio calculation (y=0) 
and with the phenomenological addendum (y=0.06) are 
shown. To quantify the accuracy of calculations, just as for 
the pairing gap in Sect. 2, we calculated the rms deviations 
of the theory from experiment, 



D 



(25) 



with obvious notation. For the lead isotopes, we find 
<5M 2/) (y=0)=0.64 MeV and(5M 2/) (y=0.06)=0.42 MeV, N p =26. 
For the tin isotopes, we get £M 2; ,(y=0)=0.68 MeV and 
<5M 2 ,,(y=0.06)=0.45 MeV, N p =2 1 . For isotones with N=82 
we obtain SM^(y=0)=0.66 MeV and 5A^(y=0.06)=0.42 
MeV, N„-19. We included into analysis only nuclei with 
the experimental error of the double difference less than 
0. 1 MeV The total disagreement for neutrons and protons 
is £M tot (y=0)=0.66 MeV and £M tot (y=0.06)=0.43 MeV. 
Thus, introducing the phenomenological addendum makes 
the agreement with data better but not so much as for the 
gap- 
One can see that the agreement is worse for lighter iso- 
topes and isotones which are close to the drip line, the 
chemical potential ju in Eq. ( l23l l being small. The reason 
is that in vicinity of drip line the concept of energy in- 
dependent mean field becomes erroneous 0211221 . If we 
limit the analysis to nuclei with fx<— 4 MeV, we obtain 
^^(y=0)=0.53MeVand^^(y=0.06)=0.27MeV,AT p =10 
for lead isotopes and <5M 2p (y=0)=0.64 MeV and 6M 2p (y-0.06) 
=0.40 MeV, N p -16 for the tin chain. For isotones N=82 

we ha ve 6M 2 „(y=0) = 0.71 MeV and 
<5M 2 „(y=0.06) = 0.39 MeV, N„=16. It is worth mention- 
ing that for the Dt quantity, Fig. [8] inclusion of the phe- 
nomenological addendum makes agreement significantly 
worse. The total disagreement for neutrons and protons is 
now <5M tot (y=0)=0.64 MeV and (5M tot (y=0.06)=0.36 MeV. 
Thus, after the inclusion of the phenomenological correc- 
tion, we obtain the agreement only a little worse than for 
the superfluid subsystem. 
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Fig. 10. Phonon corrections to the interaction amplitude f 



We observe that our calculations, which reproduce the 
value of D2 on average, do not catch the trend of the A- 
dependence of this quantity. We can suppose that it could 
be due to the fact that phonon corrections, included into 
the parameter y, are not considered explicitly. Estimations 
show that for this problem they are not so universal as in 
the gap problem as far as the diagonal element CV e ff)n m 
Eq. (l24l) for (r'(E s ))\ \ often dominates. In such a situation, 
the phonon corrections could strongly depend on the state 
H>. 

The main corrections to Eq. (1231 caused by the L-phonon 
are shown in diagrams of Fig. [10] The first one is the so- 
called induced interaction, the second one is the "ends cor- 
rection". In magic and semi-magic nuclei, the vertex g as 
a rule can be considered as a small parameter, and the so- 
called ^-approximation can be used. The ^-corrections 
are displayed in Fig. [TO] The wavy line corresponds to the 
phonon ©-function, D((l>)=2(l>l/((l> 2 -cof), where a>i is the 
excitation energy of the L-phonon, and the transferred en- 
ergy w=0 for the diagonal matrix element. The vertex g(r) 
of creating L-phonon and the frequency ojl were found 
Il23l by within the self-consistent TFFS fl24]|. Note that 
all low-lying phonons possess the normal parity n={—\) L , 
therefore we do not show it explicitly. The first diagram 
can be easily evaluated: 



is valid. However, sometimes the "resonance" cases oc- 
cur when the quantity Se\i is anomalously small. Then the 
(^-approximation does not work and higher order terms 
should be summed up. In such cases, it is natural to mul- 
tiply by y[Z\ also the ends of the first diagram in Fig. [TO] 
Then the expression for the double mass difference with 
phonon correction is as follows: 



D 2 = (D 2 + 6§D2)ZI 



(28) 



Let us calculate from this formula the corrections to 
the proton double mass difference D 2p first for the double- 
magic 208 Pb nucleus. In this nucleus the collective 3~-phonon 
with the energy W3 =2.684 MeV plays the main role. Adding 
two proton particles, we have D^iy = 0)= - 2.35 MeV, 

\l)=lh 9/2 , Z 1= 0.98, S^D 2p =-Q.U MeV and S 2/! =-2.38 
MeV. Thus, in this case the phonon correction, <5 p hL) 2p = 
L) 2 p-D 2p =-0.03 MeV is negligible. Adding two proton 
holes, we have D~ p (y = 0)= + 1.43 MeV, |l)=3ii/ 2 , ob- 
viously, we have S ^D 2p =0 due to the angular momentum 

conservation. Then, we have Z\ =0.96 and D 2p -\ .32 MeV, 
thus the phonon effect 5 p hD2 ; ,=-0.1 1 MeV is more than 
for the 1/29/2 state but also less than 10%. Another situ- 
ation takes place for the semi-magic nucleus 204 Pb. Here 
the phonon corrections appear mainly due to the low-lying 
2 + -phonon, w 2 =0.882 MeV. In this case, for the state lhg/ 2 
we have Zi=0.74, 6 ( ^D 2p =-0.55 MeV, 5 2p =-1.53 MeV 
and (5 p hD^=0.72 MeV. For the state 3si/ 2 again we have 
^D 2; ,=0, but Zi =0.82, S 2p =0.93 MeV and (5 ph D 2p =-0.47 
MeV. Thus, for both states the phonon correction is ap- 
proximately of the same magnitude as the effect of the 
phenomenological addendum in Eq. (fT3l at y=0.06. Hence 
now we must take smaller value of y. In particular, for the 
state 3si/2 the value of y=0 looks preferable. Evidently, 
more consistent approach is necessary with systematic ac- 
count for the phonon corrections with a new readjustment 
of the parameter y. 



M.) n 2((j 1 h\\Y L \\j 1 h)g u f 



to L (2ji + 1) 



(26) 



where ( \\Y[}\ ) denotes the reduced matrix element [25 1, 
and gn is the radial matrix element of g(r). 

As to the second diagram, one has the same correc- 
tions for each of the ends. They can be readily summed 
[ 1 3 1 . As the result, each of the ends should be multiplied 
by y[Z\, where Z\ is the residue of the Green function G(s) 
at the point s-ei . The corresponding correction to the dou- 
ble mass difference is equal to 



Sglh = (zt - 1)D 2 . 



(27) 



This expression contains terms of higher order in g 2 but 
such partial summation has clear physical meaning. In ad- 
dition, sometimes the difference Z\ — 1 is not small. The 
matter is that the ratio g-gnlSen, where 6e\2=£\ -E2+10L, 
is the actual dimensionless "small" parameter. For collec- 
tive L-phonons the matrix elements gn — 1 MeV, whereas 
the typical value of the denominator is 6e\2 — sf/A 1 ^ 3 , 
£f being the Fermi energy and A the mass number of the 
nucleus under consideration. Thus the condition g «c 1 



4 Conclusions 

Two aspects of the problem of evaluating double odd-even 
mass differences D 2 in semi-magic nuclei are analyzed. A 
semi-magic nucleus contains two components with differ- 
ent properties, a superfluid subsystem and a non-superfluid 
one. For the superfluid subsystem, D 2 is supposed to be 
equal to 2A, the gap A being solution of the gap equation. 
Such equalization contains an inherent inaccuracy for A of 
the order of — 0.1 -r 0.2 MeV [9] and an additional one due 
to particle number non-conservation effects 0.1 MeV 
Ifl2l . For D 2 value they should be multiplied by a factor 
of 2. In Ref. [9| the gap equation was solved with the ef- 
fective pairing interaction of the semi-microscopic model 
in which the main term found from first principles is sup- 
plemented with a small phenomenological addendum con- 
taining one phenomenological parameter y supposed to be 
the same for all medium and heavy atomic nuclei. The data 
for 34 nuclei were analyzed for isotopes of the lead and 
tin chains and isotones = 82. Neighbors of the double 
magic nuclei were excluded, as for them particle number 
non-conservation effects are especially large Ifl2l . Overall 
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disagreement of the theory with the data is -J(Jth ~^exp) 2 - 
0. 1 3 MeV for the DF3 functional flJD and =0. 14 MeV, for 
the DF3-a one 1191 . For the mass differences D2 it corre- 
sponds to accuracy a little better then 0.3 MeV. 

For the non-superfiuid subsystem, D2 is found solv- 
ing the equation for two-particle Green function of normal 
systems. This equation contains exactly the same effective 
pairing interaction as the gap equation. Results of calcu- 
lations for the lead and tin isotopes and N = 82 isotones 
with the DF3 functional and with the same value y=0.06 
of the phenomenological parameters are presented here in 
Figs. [4]-|9] It is worth to note that the problem for non- 
superfiuid nuclei, in principle, does not contain so serious 
inaccuracies as the pairing problem. The main approxima- 
tion which is made is that the two-particle Green function 
K which concerns the array ofA,A + l,A±2 nuclei is 
described in the basis of the A-nucleus. Estimation of the 
accuracy of this approximation is ~ 2 /A. However, the 
agreement turned out to be a little worse than for super- 
fluid case, (5M to t=0.36 MeV. The reason is, evidently, in 
phonon corrections which are in this case not so regular as 
in the gap equation. Estimation of these corrections shows 
that, indeed, a more consistent approach is necessary with 
systematic account for the phonon corrections, the param- 
eter y being readjusted anew. 
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